Gluon Saturation Effects in Relativistic U+U Collisions 
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We examine entropy production in relativistic U+U collisions on the basis of a Color Glass Con- 
densate (CGC) type picture as implemented in the Kharzeev-Levin-Nardi model (KLN). In this 
framework, we find that the peak entropy density produced in tip-on-tip U+U collisions is about 
30% greater than that seen in central Au+Au collisions. Although the resulting difference in the 
produced charged particle multiplicity between tip-on-tip and side-on-side collisions is smaller than 
that predicted by previous Glauber model estimates, it is still large enough to allow for experimental 
discrimination between average orientations of the uranium nuclei. We also point out that in the sat- 
uration/CGC approach the collision geometry plays a more important role than previously believed, 
and that the observed centrality dependence of the produced particle multiplicity per participant in 
Au+Au collisions can be qualitatively reproduced even without running coupling effects. 



PACS numbers: 25.75.-q, 25.75.Nq, 12.38.Mh, 12.38.Qk 

Recent studies [ij, |j, |3| have pushed for a uranium- 
uranium (U+U) collision program at RHIC to address a 
number of unanswered questions about the quark-giuon 
plasma (QGP). In particular, it is thought J., 2^ that 
the higher entropy densities produced in U+U over those 
found in Au+Au collisions would provide a definitive test 
of the seemingly ideal hydrodviiamic behavior of the el- 
liptic flow observed at RHIC [lli|. 

Previous estimates [ll y, y| relied on a Glauber model 
calculation of the entropy densities produced in U+U col- 
lisions. These results showed a substantial increase in the 
initially produced entropy density over that seen in even 
the most central Au+Au collisions. Perhaps most im- 
portantly, a clear (~ 15%) difference in the multiplicities 
produced in tip-on-tip versus side-on-side U+U collisions 
allowed the selection of average relative orientations of 
the two uranium nuclei, by placing cuts on the charged 
particle multiplicity [ij . This remains true even if exper- 
imental inefficiencies in triggering on full-overlap U+U 
collisions, using strict spectator nucleon cuts, are taken 
into account 



However, our theoretical understanding of the initial 
particle production processes during the nuclear colli- 
sion and the conditions of the produced plasma at early 
times remains limited, and it has recently been shown 
[a, l3| that this results in significant ambiguities for the 
initial source eccentricities in non-central Au+Au colli- 
sions at RHIC. Hirano and collaborators |y,l2i found that, 
if one uses parametrizations for the initial entropy pro- 
duction that are based on the Color Glass Condensate 
model J8| of gluon saturation in large nuclei at high en- 
ergies J2| , in particular the Kharzeev-Levin-Nardi (KLN) 
parametrization [lfl|, |ll|, |l3| , one obtains initial density 
profiles in the transverse plane which are flatter and have 
steeper edges than the more Gaussian-looking profiles 
[13 one obtains from the Glauber model. As a conse- 
quence, one obtains initial source eccentricities which, 
over most of the impact parameter range, are signifi- 



cantly larger in the KLN model than in the Glauber 
model. As shown in Ref . 7] , this ambiguity in the initial 
conditions has considerable impact on the interpretation 
of available elliptic flow data and the question to what 
extent the QGP created at RHIC behaves like a "perfect 
fluid". 

In the present paper we investigate to what extent 
these model ambiguities for the initial entropy produc- 
tion also affect the usefulness of a U+U collision pro- 
gram where one tries to exploit the intrinsic deformation 
of the colliding nuclei to generate deformed fireballs even 
in central collisions. Specifically, we study the influence 
of various model assumptions about the initial entropy 
production on the distribution of source eccentricities in 
U+U collisions, in order to assess how model ambiguities 
affect our ability to select for specific collision geometries 

by cutting event samples on the number of spectator nu- 

pn n n, 
cleons and on charged particle multiplicity 1, 2, 3]. 

We begin by outlining our procedure for calculating the 
initial entropy production from the Color Glass Conden- 
sate model. We use the KLN parametrization |lOl lll| , 
in the specific form implemented by Hirano and Nara in 
Ref. [13I, where the initial energy per unit rapidity and 
per unit transverse area is given by the /cT-fa-ctorization 
formula J2I] 
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with 2;i_2 — Pt exp(±y)/-\/s. The unintegrated gluon 
distribution for nucleus A is taken to be 
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where the saturation momentum Qg is determined by 
solving 

Ql{x, xj_) = -T^asiQl) xG{x, Ql) npart,A(a;±), (3) 
Of 

and similarly for nucleus B. The normalization k in 
Eq. (jSJ is fixed by the measured charged particle mul- 
tiplicity dNch/dy in central Au+Au collisions at ^/s = 
200AGeV (see below). To reproduce the results of 
Refs. and [13, we take k^ = 1.8. 

Both the energy density 0J and the local (i.e. a;^-de- 
pendent) saturation momenta of the two nuclei, Eq. Ipj)l. 
depend on the impact parameter of the collision through 
the density of wounded nucleons ripart in the transverse 
plane, computed from the standard Glauber model for- 
mula (see, e.g., Ref. ll3J). For nucleus A it is given by 
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while for nucleus B the nuclear thickness functions Ta 
and Tb should be interchanged. The latter are computed 
from a Woods-Saxon distribution for the density of the 
uranium nucleus, with radius R{9) — (6.8 fm)(0.91 -f- 
0.26 cos^ 6*) (where is the polar angle with the nuclear 
symmetry axis) and surface thickness ^ = 0.54 fm. 

The integrated gluon distribution, xG{x, fc^), in Eq. (0 
is given by 



xG(x, k 
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with K = 0.7 and A = 0.2 [lj| . The running coupling con- 
stant tts is computed from the one-loop expression 
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with /3o = 11— §"■/ = 9. To render the calculations in- 
frared safe, we limit the growth of the running coupling 
constant at small k^ by cutting it off at as{k^)=0.5 
In all expressions above we use A = Aqcd = 0.2 GeV 




From equation (^ we obtain the energy density at 
thermalization time To (for which we use tq = 0.6 fm/c) as 
e{x±, To) = dE/{To d^x_L drj) at ?7 = y = 0. We then trans- 
late this energy density e into an entropy density s, by 
assuming that it thermalizes and using the relationship 
s = 8.67e'^/^ (where e must be inserted in GeV/fm'^) for 
a thermalized ideal gas of quarks and gluons. By adjust- 
ing K, this entropy density is then normalized such that 
for impact parameter 6 = its integral over the trans- 
verse plane yields the same result as the Glauber model 
in reference [l| (c/. Figure^. 




y(fm) 



FIG. 1: Transverse entropy density profiles in the in-plane 
direction (left) compared with those in the out-of-plane direc- 
tion (right) for central Au-|-Au (top panels), tip-on-tip U+U 
(center panels) and side-on-side U-l-U collisions (bottom pan- 
els) . Shown are results from the Glauber model (dotted) and 
from the KLN model with three different pr cutoffs (see text): 
Pcut = 3 GeV/c (dashed gray), pcut = 2 Q™"" (solid black), and 
Pcut = oo (solid gray). 



We then assume that the charged particle multiplicity, 
dNch/drj, is proportional to the total entropy produced 
in the plane, dNch/drj = T J d'^r±s{r±), and adjust the 
proportionality constant F = 0.088 to agree with cen- 
tral Au-|-Au data obtained by the PHOBOS collabora- 
tion [lj|. The resulting fit is shown in Figure |2| These 
parameters are then assumed to also describe U-l-U col- 
lisions at the same collision energy. 

In the work of Hirano and Nara, the integral over px 
in Eq. (Q) is cut off at a fixed value of Pcut = 3 GeV, to 
limit the contribution from high-p^ gluons which ther- 
malize incompletely ^2i ■ As shown in Fig. |21 the result- 
ing centrality dependence of the charged particle mul- 
tiplicity in 200 A GeV Au-|-Au collisions reproduces the 
PHOBOS data [1J| very well. However, the use of a fixed 
Pt cutoff in Eq. (Q, independent of impact parameter, 
raises questions: The momentum dependence of the un- 
integrated gluon distributions <^a,b is controlled by their 
respective saturation momenta in Eq. ||2Il which do de- 
pend on both x± and b; by implementing a fixed cutoff in 
Pt, one doesn't allow the px integral to fully explore this 
6-dependence. As we discuss now (and show in Figure 
I^J, this leads to a significant distortion of the centrality 
dependence of charged hadron production. 

At any given point x±_ in the overlap region of the 
two nuclei, the result of equation ^ is controlled by 
two saturation momenta, Qs^min and Qs^max, correspond- 
ing to the smaller and larger of the nuclear thicknesses 
of the two colliding nuclei at this point. If we ignore 
the running of as, the integrand of equation ^ is es- 
sentially constant for pt <Qs,min, decreases like l/_p^ 

for Qs,min <PT< Qs.max and fikc 1/p^ for Pt > Qs.max- 

With a fixed px cutoff we cut off more of the l/p^ tail 
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FIG. 2: Charged particle multiplicity per participant pair as 
a function of centrality for Au+Au. The inset shows cor- 
responding predictions for full-overlap U-fU collisions. The 
data are from PHOBOS at y5 = 200A;GeV [H]. Shown 
are results from the Glauber model (dotted) and from the 
KLN model different py cutoffs: pcut = 3GeV/c (dashed), 
Pcut = 2 Q™'"' (black), and pcut = oo (gray). 



in central than in more peripheral collisions, due to the 
fact that Qs^max is larger in the former case. The net 
effect of this is to suppress entropy production in cen- 
tral relative to peripheral collisions. Indeed, the central- 
ity dependence of dNch_l dy becomes noticeably steeper 
if this pt cutoff is removed (solid gray line in Figure 
121). The same effect can be achieved by allowing the 
Pt cutoff to scale with the (centrality dependent) effec- 
tive saturation momentum. We tested this by setting 
Pcut = 2Q™'''^(6) where we chose (somewhat arbitrarily) 
Qfa'^(6) = Q,(j;=0.0I,a;±=O;b) (solid black line in Fig- 
ure [2Jl. 

Even if for Au+Au collisions the effect of removing the 
Pt cutoff in Eq. ^ is clearly noticeable in the centrality 
dependence of dNch/dy, it remains quantitatively small 
(~ 15% in peripheral collisions). For full-overlap U+U 
collisions, however, the different ways of cutting off or not 
cutting off the pt integral produce qualitatively different 
results. The inset in Figure |21 shows that for a fixed pt 
cutoff of 3 GeV/c (dashed black line) , tip-on-tip collisions 
produce lower multiplicities than side-on-side collisions, 
in stark contradiction with the Glauber model prediction 
(dotted line). The qualitative tendency of the Glauber 
model toproduce higher multiplicities in tip-on-tip colli- 
sions 111, y, y| (due to the larger binary collision compo- 
nent in this configuration) is recovered in the KLN sat- 
uration model if the pt cutoff is either entirely removed 
(solid gray line) or assumed to scale with the saturation 
momentum (solid black line). In either case, however, 
the quantitative difference between tip-on-tip and side- 
on-side multiplicities is somewhat smaller for the KLN 




FIG. 3: Eccentricity versus impact parameter for Au+Au col- 
lisions at y/s = 200 GeV for the Glauber model (dashed) and 
the KLN model with pcut = 2Qs (black), pcut = cx3 (gray) and 
Pcut = 3 GeV/c (dashed gray) . 



model than for the Glauber parametrization. 

Returning to the initial entropy density profiles shown 
in Figure ^ we can now see how these differences arise: 
Generically, the profiles from the saturation model are 
flatter near the center and steeper near the edge than 
those from the Glauber model (which look more Gaus- 
sian) . In tip-on-tip collisions (center panels) the Glauber 
model shows the strongest growth of entropy density in 
the middle of the overlap region, due to a strong increase 
of the binary collision component. This growth is tem- 
pered in the KLN saturation model, due to the slower 
(logarithmic) increase of the saturation momentum Qg 
as one approaches the center of the overlap region. If one 
implements a flxed pt cutoff, the pt integral in Eq. (Q 
can not even fully explore this logarithmic growth of the 
saturation momentum, and the entropy density increases 
even more slowly towards the center. In fact, the differ- 
ence in entropy density between tip-on-tip and side-on- 
side collisions then becomes so small that it can no longer 
compensate for the smaller overlap area in tip-on-tip col- 
lisions, leading to the aforementioned counterintuitive re- 
sult of smaller total multiplicities in tip-on-tip compared 
to side-on-side collisions (dashed gray line in the inset of 
Fig. H. 

As already pointed out in Refs. 0,13 > the more plateau- 
like entropy distributions from the KLN model yield 
larger eccentricities €x — ('i.^'il than the Glauber model 
profiles, due to their bigger weights at larger distance 
from the fireball center. In Figure 13 we show the eccen- 
tricity of the initial entropy distribution as a function of 
impact parameter b for Au+Au collisions. We confirm 
the observation made in [Tf that the KLN eccentricities 
exceed the Glauber ones by about 25-30% over most of 
the impact parameter range, and further show that they 
are insensitive to the choice of the pt cutoff in Eq. JQ. 
Since the value of the initial source eccentricity controls 
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FIG. 4: Multiplicity distribution for full-overlap U+U col- 
lisions. Shown are results from the Glauber model (dotted) 
and the KLN model with pcut ~ 2 QJ"*" (dashed) and pcut ~ oo 
(solid black). Also shown is the result of cutting the full mul- 
tiplicity distribution given in Figure |H] on the 0.5% of events 
with the lowest number of spectators. 



FIG. 5: Eccentricity distributions resulting from cutting the 
multiplicity distributions of Fig. |l|on the 10% of events with 
the lowest (black) and highest (gray) multiplicities. Shown 
are results for the Glauber model (dotted) and the KLN model 
with Pcut = 2 Q™ ^^ (dashed) and Pcut ~ oo (solid) . 



the finally observed elliptic flow, Fig. |31 reassures us that 
the conclusions of Ref. |3 as to the interpretation of the 
measured V2 values are not affected by model uncertain- 
ties related to the choice of pr cutoff. 

In [l|, y two of us advocated using the multiplicity to 
select for particular orientations of the colliding uranium 
nuclei. In the Glauber model framework, we found that 
of all possible orientations of full-overlap U-l-U collisions, 
the highest multiplicities were produced in tip-on-tip col- 
lisions and the lowest multiplicities were formed in those 
collisions between nuclei in the side-on-side configuration. 
The difference in mean multiplicity between these ex- 
treme cases was ~ 14% [l|. This provided a "handle" to 
select for average orientations of the nuclei: By selecting 
events with high multiplicity, one effectively triggers on 
low source eccentricity, and vice versa. The 14% differ- 
ence in multiplicities between tip-on-tip and side-on-side 
collisions proved large enough that such a discrimina- 
tion between collision configurations with large and small 
source eccentricities remained feasible even after exper- 
imental inefficiencies in triggering on full-overlap U+U 
collisions (by selecting events with a small number of 
spectator nucleons) were taken into account ;2'| . 

The inset in Fig. El shows that with CGC/KLN initial 
conditions the range of multiplicities produced by the 
different orientations of full-overlap collisions is some- 
what smaller than in the Glauber model calculations. 
It may therefore be prudent to reanalyze the possible 
degradation of our ability to select event classes of well- 
defined source eccentricity resulting from trigger ineffi- 
ciencies and event- by- event multiplicity fiuctuations. 

We first study the effect of event-by-event multiplicity 
fiuctuations, assuming that full-overlap collisions can be 
selected with 100% efficiency. As before [l| we adopt as 
probability distribution for the event-by-event multiplic- 



ity the expression [l(]| 
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where n is the measured multiplicity in a given event and 
n(^) is the average multiplicity for a given orientation $ 
of the nuclei, now computed from the KLN model us- 
ing the procedure described above. For the parameter a 
controlling the width of the fluctuations we again take 
a = 0.6 as this has been shown to produce good agree- 
ment with existing Au-|-Au data (lOl]. 

Integrating Eq. ((2j) over $ yields the multiplicity distri- 
bution for full-overlap U-l-U collisions, shown in Fig. 01 
As expected, the distributions obtained from the KLN 
model are considerably narrower than the one resulting 
from the Glauber model. Furthermore, the shapes of the 
distributions arc quite different, a result of the individual 
Gaussian distributions from tip-on-tip and side-on-side 
collisions being squeezed on top of one another, due to 
the smaller difference in mean multiplicity between them. 

After placing cuts on the multiplicity distributions of 
Fig. ^ we can convert the multiplicity distribution into 
a distribution of eccentricities according to 
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where no and tii represent the lower and upper limits 
of the multiplicity cuts, respectively. The resulting ec- 
centricity distributions are shown in Fig. O for cuts on 
events with the 10% lowest and highest multiplicities, re- 
spectively. Not unexpectedly, the distributions resulting 
from the KLN model are slightly broader than those from 
the Glauber model: Due to the smaller range of multi- 
plicities between the $ = and <& = § configurations, the 
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FIG. 6: Multiplicity distribution for U+U collisions with arbi- 
trary impact parameters and relative orientations. The gra 
curve corresponds to the Glauber model calculation of Ref. 
while the result from the KLN model is given by the solid 
curve. 



eccentricities corresponding to given multiplicity slice are 
less sharply defined. This also results in a slightly higher 
(lower) average eccentricity of high (low) multiplicity 
samples: For the high multiplicity cut, the KLN model 
(with the pt cutoff choices listed in Fig. \^ gives an av- 
erage eccentricity (cx) =0.03-0.04 (vs. (ex) =0.014 from 
the Glauber model), whereas for the low multiplicity cut 
the KLN model gives (e^) =0.22-0.23 (vs. (£3;)= 0.24 
from the Glauber model). We conclude that the differ- 
ences between the two classes of models are small, and 
the suggested event selection scheme remains feasible in 
this new framework. [It is worth noting that we have 
computed these same eccentricity distributions for en- 
ergies that will be available at the LHC and found no 
significant deviations from the results presented here.] 

We now proceed to also discuss the effects of trigger 
inefficiencies in the selection of full-overlap collisions on 
the eccentricity selection. We repeat the calculations re- 
ported in Ref. |2| for the KLN model, restricting our at- 
tention to the case without pr cutoff {pcut — *■ 00). Figure 
El gives the resulting complete multiplicity distribution 
for U-|-U collisions of any impact parameter and relative 
orientation, and compares it with the previously studied 
Glauber model initialization. The most notable feature 
is the reduced width of the high-multiplicity tail in the 
KLN model, due to the smaller spread of multiplicities in 
6 = 0, full-overlap collision events seen in Fig. 01 The lat- 
ter can be selected from this distribution by placing very 
strict cuts on the number of spectator nucleons. The 
gray histogram in Fig. ^ shows the multiplicity distri- 
bution for an event class that was obtained by selecting 
from the KLN distribution in Fig.Othe 0.5% events with 
the lowest spectator counts. If we cut this multiplicity 
distribution once more for high (low) multiplicities, we 




FIG. 7: Eccentricity distributions resulting from cutting the 
gray histogram in Figure|4lon the 5% of events with the lowest 
(black) and highest (gray) multiplicities. The corresponding 
average eccentricities are 14% (black) and 4% (gray). 



can again select events with small (large) average eccen- 
tricity. 

The resulting eccentricity distributions are shown in 
Fig. [7| The black (gray) histogram corresponds to the 
subclass of events obtained by selecting from the gray 
histogram in Fig. 01 the 5% lowest (highest) multiplicity 
events. While the eccentricities are rather widely dis- 
tributed within these two event subclasses, their aver- 
age eccentricities are very different: (£3;)= 0.04 for the 
high-multiplicity subclass vs. {ex)=0-lA for the low- 
multiplicity selection. (For the Glauber model the cor- 
responding average eccentricities were 0.02 and 0.18, re- 
spectively, and actual values distributed somewhat more 
narrowly about these averages p|.) 

We conclude that, in spite of the documented quan- 
titative differences between the initial entropy distribu- 
tions calculated from the KLN and Glauber models, the 
proposed event selection scheme for U-l-U collisions con- 
tinues to allow isolating event classes with widely differ- 
ent source eccentricities, at energy and entropy densities 
that exceed those reachable in Au-|-Au collisions by a 
wide margin. Detailed comparisons of elliptic flow mea- 
surements in U+U collisions with ideal fluid dynamical 
calculations will, however, be affected by these model un- 
certainties. 

Before closing, we would like to point out the impor- 
tant role played by the collision geometry in determining 
from the saturation model the centrality dependence of 
particle multiplicity per participant. For simplicity, let 
us assume that the multiplicity of produced gluons at a 
given location x± in transverse plane is given by 



dN 
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mm{Ql^{x^),Q''^^{x^)} 



(9) 



Here Qli{x±) and (5^2 (^i) ^^^ ^^e saturation scales in 
the two colliding nuclei at the given transverse position 




400 



FIG. 8: Particle multiplicity per participant as a function of 
centrality, for the simple CGC-motivated geometrical model 
Eq. Q. The proportionality constant in Eq. Q is chosen to 
fit the most central data point from |14| (see text for details). 



x^. Equation 1^1 captures with logarithmic accuracy the 
main quahtative behavior of integrated particle multi- 
phcity in the saturation/CGC approach. It can be ob- 
tained from Eq. IpJ if one neglects logarithms of satura- 
tion scales as slowly- varying corrections to the power-law 
scaling of Eq. © . 

For illustration, let us consider a collision of two spher- 
ical nuclei (e.g. a Au-|-Au collision). If the radius of the 
nuclei is R and they collide at impact parameter 6, the 
saturation scales are (assuming box-like nuclear density 
distributions for simplicity) 



QI^{X^) ex .^i?2 - {x^ + \bY, 

Ql^ix^) ex ^R^-{x^-\bY, 



(10) 



with the same proportionality coefficient for both of 
them. Remembering that the square of the satura- 
tion scale is proportional to the density of participants 
npia:t{xi_), we deduce that the initial gluon (and thus also 
the final charged hadron) multiplicity per participant is 
proportional to 
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with the total number of participant nuclcons 
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(X d^x^{Ql,{x^) + Ql^{x^)). (12) 



The integrals in Eqs. (|ll|l and (|12(l with the saturation 
scales Hl()|l can be done numerically, yielding the plot 
of particle multiplicity per participant as a function of 
Npart shown in Fig. |H1 The normalization of both axes 



was adjusted to fit the most central data point in Fig. |5] 
The curve in Fig. |S1 is seen to qualitatively reproduce 
the monotonous rise of (\ / N-pait)dN c\i/ dy with A^part ob- 
served at RHIC nil (see Fig. EJ. We stress that the 
positive slope of the curve in Fig. |S| arises entirely from 
the geometry of the colliding nuclei, as implemented in 
the simple CGC-inspired model Eq. ^ for particle pro- 
duction. It is likely that a similar geometrical mechanism 
is responsible for the positive slope of the gluon multi- 
plicity per participant obtained in numerical simulations 
of classical gluon fields performed in [I3. Although of- 
ten stated otherwise, running coupling and/or other log- 
arithmic effects appear not to be necessary to obtain this 
positive slope, but they seem to be required to turn the 
qualitative agreement of Fig. |H1 with the data in Fig. [3 
into a quantitative one. Still, we are amazed how well 
the data shown in Fig.|2]can be qualitatively understood 
within the simple geometrical saturation model of gluon 
production given in Eq. JHJ. 

Although a number of different questions can be stud- 
ied with U+U coUisions [l|, |la, 113, llal , perhaps the most 
persuasive argument in favor of a U-l-U program at RHIC 
is the potential |]| to test the ideal hydrodynamic behav- 
ior of the elliptic flow at higher entropy densities and with 
larger fireballs than are available in even the most cen- 
tral Au-|-Au collisions. In the Glauber model calculations 
presented in [ij, |2| we found that the maximum entropy 
density available in tip-on-tip U-l-U collisions was approx- 
imately 40% greater than that seen in central Au+Au. 
Examination of Figure ^ leads to a similar conclusion 
for the KLN model results, with increases of 35% and 
29% for the cases with variable and infinite px cutoff, 
respectively. Since the initial maximum entropy density 
scales roughly (although not exactly [3|) with the observ- 
able {l/{S)){dNch/dy), where (S) is the mean transverse 
overlap area of the two colliding nuclei, this increase in 
initial entropy density provides an important lever arm 
for testing the approach towards ideal fluid dynamical 
behavior at RHIC experimentally. The additional pa- 
rameter space opened up by U-l-U collisions at RHIC 
could prove even more important if Pb+Pb collisions at 
the Large Hadron Collider LHC (which are expected to 
produce much higher initial entropy densities and tem- 
peratures) turned out to produce a more weakly coupled 
quark-gluon plasma than RHIC collisions, and the QGP 
created at RHIC thus were to yield the best possible ap- 
proximation of a "perfect fluid" . In this case the U-fU 
collision program could cover a key window in parame- 
ter space, but we also note that fully exploiting it would 
depend on our ability to gain sufficient control over the 
presently remaining model uncertainties for initial parti- 
cle production, as studied here. 
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